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Charged particles in static electric and magnetic fields have Landau levels and tunneling states 
from the vacuum. Using the instanton method of Phys. Rev. D 65, 105002 (2002), we obtain the 
formulae for the pair-production rate in spinor and scalar QED, which sum over all Landau levels 
and recover exactly the well-known results. The pair-production rates are calculated for an electric 
field of finite extent, and for the Sauter potential, both with a constant magnetic field also present, 
and are shown to have finite-size effects. 
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I. INTRODUCTION 



Vacuum polarization and pair production are two physically important phenomena of quantum electrodynamics 
(QED) in strong electromagnetic fields pJQ- The one-loop effective action in constant electric and magnetic fields 
has nonlinear contributions to the classical action and, thereby, leads to the nonlinear Maxwell equations. Another 
interesting phenomenon is the pair production due to vacuum instability in the presence of electric fields near and 
above the critical electric field strength E c = m 2 c 3 /eh (1.3 x 10 16 V/cm) 0,0 (for references, see, e.g., @> lib- 
Strong QED has many physical applications In particular, electromagnetic fields of some neutron stars and 
black holes [6( and high- intensity laser fields [j| above the critical strength E c /c = B c = m 2 c 2 /eh (4.4 x 10 13 G) have 
revived recent interest and applications of strong QED. The one- loop QED effective action is known exactly in the 
background of constant electric and magnetic fields Q|. In the case in which E 2 — B 2 and E • B are not both zero, 
one can go to a frame in which E and B are parallel with magnitude E and B, and obtain the imaginary part of the 
one-loop effective action per four- volume for spinor QED given by 0, 0> H @ 



• M = W W _ (lE)(<lB) ^ 1 ^ (rmB 

■"spinor sp 

and that for scalar QED by [HH [HG3 

„ T „m ID (qE)(qB) ^ (— 1)" +1 , /mrB\ f nnm 2 \ 
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The pair-production rates per volume per time themselves are the first terms in each of these series [T^ , 
and 

,2 



Kruglov derived a general formula for the more general case of particles of arbitrary spin s with ADM (anomalous 
magnetic moment) and EDM (electric dipole moment) p^ |. 
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In this paper we apply the recently-introduced instanton method |3( to find the pair-production rate in constant 
electric and magnetic fields in spinor and scalar QED, agreeing with the results above. The idea of the instanton 
method, the elaboration of the role of tunneling in pair production 0], is that for static electric fields, fermions and 
bosons have tunneling states from the vacuum in the Coulomb (space-dependent) gauge, and their pair-production 
rates arc determined by these instanton actions for tunneling. In the presence of static magnetic fields, charged 
fermions and bosons have discrete spectrum of Landau levels. Taking into account both Landau levels and instanton 
actions for tunneling, we obtain equivalent formulae for pair-production rates in spinor and scalar QED. These formulae 
are given as sums over all Landau levels, and we recover exactly the well-known results. Further, applying the method, 
we estimate the pair-production rates by an electric field with a finite extent and a constant magnetic field, which 
have additional factors determined by the potential difference across the boundary beside the pair-production rate 
for a constant electromagnetic field. 

The organization of this paper is as follows. In Sec. II we generalize the instanton method to the case of constant 
electric and magnetic fields and obtain the pair-production rate as the sum of Landau levels. We also compare 
our instanton method with the instanton method in Euclidean time. In Sec. Ill we apply the instanton method 
to inhomogeneous electric fields: a constant electric field confined to a finite region and one from the Sauter-type 
potential. 

II. INSTANTON METHOD FOR PARTICLES WITH SPIN IN ELECTROMAGNETIC FIELDS 

With the spin properly taken into account, the component equation for minimally-coupled particles with spin in a 
constant electric field takes the form (in units with % = c = 1 and with metric signature (+,—,—,—)) 



where q (q > 0), m and a are the charge, mass and spin of the particles. We shall show the physical implication of 
the imaginary part from the spin effect. After mode-decomposition, 



rf{% + iqA^){d v + iqA v ) + m 2 + 2iaqE $ CT = 0, 



(5) 



-d 2 -{uj + qEzf + m 2 + k 2 ± + 2iaqE $ CTwki (z) = 0. 

The wave function © has a solution in terms of the complex parabolic cylindrical function E(x,y) 

$^k ± (z) = CE(a s ± x ,£), 



(6) 



(7) 



where 




(8) 



In the two asymptotic regions the wave function becomes 



A^ ± (0 ~ B^ ± (£), (£ « -2|5 sk J 1 / 2 ) 
CV™k.(0> (C»2|a sk J 1/2 ), 



(9) 



where 




(10) 



and 



A = iC(l + e 2 ™° k ±) 1/2 1 B = ~iCe 



: _L 



(11) 




(12) 



The reflection probability is given by 



3 



where S is the instanton action 



o 



This is with the boundary condition that the current flux on the right hand side is outward. However, then the group 
velocity there is inward (to the left). If we instead impose the boundary condition from causality that signals are 
outgoing on the right, then on the left hand side the outgoing flux becomes \A/B\ 2 = 1 + e~ 2S times the incoming 
flux, an amplification by the Klein-paradox. Then by the results of Nikishov |l3j, the pair-production rate is just this 
amplification factor minus 1, or e~ 2S . It is interesting that the instanton approximation with positive action S (exact 
for a uniform field) never gives more than 1 expected pair per mode, though there is no such restriction for a general 
electric field. 

On the other hand, for fermions {a = 1/2,3/2, •••) there is the Klein-paradox, and the flux conservation now 
becomes 

\A\ 2 + \C\ 2 = \B\ 2 . (15) 

This is also with outgoing flux but incoming signals on the right. Replacing this by outgoing signals on the right makes 
the outgoing flux on the left become |A/i?| 2 = 1 — er 2S times the incident flux on the left. (There is no amplification 
factor larger than one in the fermion case with causality imposed, because of the Pauli exclusion principle, as Feynman 
once explained to one of us (DNP) while drawing diagrams and saying, "I'm supposed to be good at these diagrams.") 
Then by Nikishov's results |l3|. the pair-production rate is just one minus this reflection coefficient, again e~ 2S . This 
result confirms the use of the instanton method in this paper and Q . 
In the instanton method 0, we use the Klein-Gordon equation 

rf^d,, + iqA, t ){d v + iqA u ) + to 2 ] $(i,x) = 0. (16) 

The significant difference between Eqs. iJSJ and l|I6|l is the imaginary constant term. We illustrate the instanton 
method by first considering the case of a pure electric field along the z direction. In the Coulomb gauge, the 4- 
potential is given by 

A fi = (-Ez, 0,0,0). (17) 
The component field of the Klein-Gordon equation, 

[(dt - iqEzf -d 2 x -d 2 -d 2 + m 2 ]$ = 0, (18) 

has a solution of the form 

$ = e i ^- x - L - ut) $ uk± (z), (19) 

where kj^ and x_l denote the momentum and the vector perpendicular to the electric field. Then the above equation 
becomes 

[-0 2 -{u + qEzf + m 2 + ki]$ wkx (z) = 0. (20) 

Note that Eq. I|2U|I describes the wave equation for an inverted harmonic potential with a negative energy. Therefore, 
one has tunneling states under the potential barrier, whose instanton actions are given by |3J 

2S kx =^Q^( 2 ) = ^!^l, (21) 

where Q(z) = m 2 + — (u> + qEz) 2 . The tunneling state and instanton action can also be understood in classical 
theory. The mass of the charged particle is invariant: 

m 2 = V^iPn - qAJfa, - qA v ) = (p + qEz) 2 - p 2 z - p\. (22) 

The time component and the transverse component of 4-momentum are constants of motion, say, po = lo and p± = kj_ . 
In the Euclidean spacetime t — it, the z component becomes p z — ipf, where pf is the Euclidean momentum. Hence 
Eq. 11221) becomes a harmonic oscillator in the reduced phase space (z,p z ), 

^ + (f! (z + i£_)^ + ki) . (23) 



4 



The action of the oscillator, given by the energy divided by the frequency (qE)/(2n), is the instanton action 

5^=2^=7^^. (24) 

The instanton action 1)2411 can also be obtained in the time-dependent gauge = (0, 0, 0, — Et). Using the action 
S = J L(t,x), Popov jl(]| obtained the same instanton action in the Coulomb gauge, and Affleck et al. a l so 
obtained the same result in the mixed-gauge A^ = —F^Xv/2. Therefore, the instanton action does not depend on 
the choice of gauge, which is just a matter of technical simplicity. 

In a pure magnetic field B = Be z along the z-direction, A = (0, Bx,0), and the component field equation of the 
Klein-Gordon or Dirac equation has a solution of the form 

$ = e i(M+*^-«<)$ ff± ^ fea ( x ) 5 (25) 

which leads to the mode equation 

-dl + {qBx + k y f +m 2 + k 2 z - 2qBa± - w 2 ] ®„ ± u>k v k z (x) = 0. (26) 

Here a± is the spin projection: a = for the scalar and a± = ±1/2 for the Dirac spinor. The mode equation l|2t)|) 
has bound states given by harmonic wave functions with the energy spectrum |18| 

uj 2 = m 2 + k 2 z +qB(2j + l-2a ± ), (j = 0,l,---). (27) 

The discrete spectrum due to the magnetic field is the Landau levels for charged particles. Note that all the Landau 
levels are non-degenerate for the scalar particles, whereas all the states of the Dirac spinor are doubly degenerate, 
j, a + = 1/2 and j — 1, cr_ = —1/2, except for the unique lowest Landau level, j = 0, cr+ = 1/2. 

Now we apply the formulae [3( for pair production based on the instanton calculation to the static uniform electric 
and magnetic fields. For the electric and magnetic fields parallel to each other along the z-direction, the 4-potential 
is given by 

Afj, = (—Ez, 0, Bx, 0). (28) 

The component field equation, 

[{d t - iqEzf -d 2 ~ (d y + iqBx) 2 - d 2 z + m 2 - 2qBa ± ]$ rT± = 0, (29) 
has the solution of the form 

® a± =e l ^y-^^ ±ujky (x,z). (30) 

Then the above equation becomes 

[{-dl + (qBx + ky) 2 } + {-8 2 -{uj + qEz) 2 } + m 2 - 2qBa ± ]<$> a±Ljky {x 1 z) = 0. (31) 
The first parenthesis has the harmonic wave functions as cigenfunctions, so the remaining equation becomes 

[-dl ~{lu + qEz) 2 + m 2 + qB(2j + 1 - 2a±)]$ (7±ukyj (z) = 0, {j = 0, 1, • • • , ). (32) 

Therefore, the charged particles, exactly described by the inverted harmonic potential, have, as the instanton action 
for tunneling states, 



2Sa±j = 7T - ' H ±; - (33) 



m 2 + qB(2j + 1 - 2a±) 
qE 

The main result of Ref. Q , when converted from the imaginary part of the effective action to the pair-production, 
is that the expected number of pairs produced per mode in a static electric field is given by 

N± = ±(e ±w - 1) = e- 2S , (34) 

where the upper (lower) sign is for bosons (fermions) and S is the instanton action for the corresponding mode, here 
given in Eq. (|33|l . We then obtain the pair-production rate [l9| for scalar particles as 

to^l^^ K'^"' 4 " ), (35) 
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and for fermions as 

«r(i) , P R s (ggKgg) ^ v-^ / m 2 + ^(2j + l-2a ± ) A 

V ' J=0 C r ± =±l/2 ^ 

Here (gE)/(2n) is the (number) density of states (per momentum) from the w-integration, and (qB)/(2ir) is that 
available for each Landau level. When we sum the geometric series of Eqs. H35fl an d H36|) . we readily get the standard 
results of Eqs. © and respectively. 

It should be noted that the pair-production rate does not depend on the renormalization scheme, since all divergence 
and renormalizability is contained in Re£ 0, IT^ | • 



III. INHOMOGENEOUS FIELD 



We now extend the analysis to inhomogeneous electric fields together with a constant magnetic field. As inho- 
mogeneous electric fields we consider a localized constant electric field in the region \z\ < L and an electric field 
obtained from the Sauter potential Aq(z) = — EL tanh(z/L) for qEL ^> m and mL ^> 1. In both cases, the electric 
field extends effectively a distance 2L in the z-direction. For an electric field localized in the z-direction, pairs are 
produced only when u) — qAo(+oo) > m and uj — qAo(-oo) < —m. Also the instanton actions exist when the Landau 
levels are limited to 

qB(2j max + 1) = mm{(w - qA a {+oo)) 2 - to 2 , (w - qA a {~oo) f - to 2 }. (37) 

The pair-production rate per area and per time (using an overbar to distinguish from rates per volume per time) now 
takes the form, for scalar QED, 

*£L = ^W ^e"^, (38) 

{2ny J qAo (+oo)+ m ~ 



and for spinor QED 



(!) (qB) r«M-oo)-m _ 

Spinor = 7o~^j dLJ l^ 1^ e ±3 l ( 39 ) 

\W) J q A (+oo)+m j= 0(7± =±l/2 



where (qB) /(2tt) is the number of state for Landau levels and another factor 1/(2tt) is from the ui integration. The 
to integration yields the potential energy difference 

V = qAo(-oo)-qA a {+oo) = 2qEL, (40) 

where L would extend to infinity for the homogeneous electric field. Here the instanton actions are determined by 



^ = 5>S' (4i) 



A:=0 



where the dominant (0-loop or classical) contribution comes from 

Sjfi = / + dz^/m? + qB(2j + 1 - 2a±) - (w - qA [z)f, (42) 

where z± are turning points of the integral and a = for scalar and o± = ±1/2 for spinor QED. 

In the first case the electric field is confined to a finite region of length 2L. The potential is given by 

Af, = (-Ez, 0, Bx, 0), \z\ < L. (43) 

The instanton actions are given by 

2S 3 = ^ + q B(2j + l) 
qE 
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for j smaller or equal to the highest Landau level 



Jr. 



2qB 



(f-M)' 



(45) 



Finally, after the lu integration, we obtain the pair-production rate per area and per time for scalar QED 

,2 . 



AT, (1) 



scalar 



___csch(— jexp(- 



x(2L) 



2m 
~V~ 



2m 
~V 



exp 



qE 

irB 
~~E~ 



y/{V/2m)*-l 



dx- 



xe i B 



(46) 



and that for spinor QED 



<1) 

spinor 



(2tt)5 
x(2L) 



(^) exp (- 



i - 



2m 2m ( nt> 



qE 



(47) 



We note that the first factor is the pair-production rate per volume for the constant electric field in Sec. II, and the 
localization of the electric field gives the factor on the second line, which is roughly 2L for V ^S> 2m. In the large L 
limit, after dividing by 2L to convert from a rate per area to a rate per volume, we recover Eqs. and fljl. since 
m/V -> 0. 

Next, we turn to the slowly varying electric field E(z) = Esech 2 (z/L) given by the Sauter potential 

(48) 



A (z) = -ELtanhy—y 

for qEL ^> m and mL 3> 1. The highest Landau level is the same as Eq. 145|) . However, the instanton actions are 
given by 



25', 



qE 



ttB[3 
E 



{2j + 1 - 2cr±) + O 



V 2 ' V 4 ) ' 



(49) 



where 



9 2 
9 1 

lu m 



(50) 



The instanton actions are also obtained from Eq. (47) of Ref. |j by replacing by qB(2j + 1 — 2a±). Then the 
pair-production rate per area and per time for scalar QED takes the form 



scalar 



(gg) 

2(2tt) 2 



c?w csch ^ ex p( 



7rm 2 a 



'-»K-^)-{-i(^-H) 2 - a )} 



(51) 



The spinor case is obtained by replacing csch(ir B f3 / E) by 2coth(7ri3/3/i?), the 2 being for two spins. In the L = oo 
limit, with J duj = V = 2qEL = oo, we recover the standard result © and (@J when we divide by 2L to get the rate 
per volume instead of the rate per area. In fact, we can see that the final factor of the right hand side of Eq. I|51|) . 
the part inside the square brackets, is nearly one over most of the integral if qEL 2 ^> 1, so that that factor can be 
dropped. 

It is interesting to compare the pair-production rates per area given by the instanton method using just the 
classical action with the rate per area given by integrating the uniform-field rate per volume over the extent of the 
inhomogeneous field: 



AT U = J dz U{E{z),B 1 



m). 



(52) 



7 



Here the subscript u is for "uniform field" . 

For simplicity, let us take the case B = 0, so the uniform- field rate per volume is the same for bosons and for 
fermions (per spin state) and in the Sauter electric field E(z) = Esech 2 (z / L) would be 

N= (^! S e C h 4 (i)expf-^cosh 2 f^l. (53) 



(2tt) 3 \L 

If we let 



qE \L 



(54) 

(where this E is the maximum value of the Sauter electric field), and if we let x = sinh(z/L), we get 



dx 



(2n) 3 ' J_oc(l + ^) 5 / 
Analogously, if in Eq. 1|51[1 we set x — 2u/V and 

2m to 



: e-*r. (55) 



< 1, (56) 



V qEL 

then assuming 5e 2 — 1/ '(nqEL 2 ) < 1 so that the final factor in Eq. (|51|l may be dropped, that equation gives 

dx 



_ (qE) 2 L _i f 1 



re 



"T" 



(57) 



Here the subscript i stands for "instanton method" , here using just the 0-loop or classical action. 
If we now have a maximum electric field E far below the critical value to 2 /q, so <5 <C 1, then 



77 u « e^ATi « Vi-L—y/rie-i = v ; ; e ?. 58 

(Z7T) J (2^)^771 



< 1, (59) 



The two estimates for the rate agree if 

6 2 7TTO 4 7TTO 4 27TTO 4 L 

45 ~~ 4(qE) 3 L 2 ~ qEV 2 ~ V 3 

but they do not agree if this quantity is not small. 

On the other hand, if the maximum electric field E is far above the critical value m 2 /q, so 5 ^> 1, then necessarily 
e 2 /S <C 1, and we get 

A7 U « ftf, « te^yfiS. (60) 

Thus in this case the uniform-field approximation gives roughly 85% of the 0-loop instanton-method approximation 
using only the leading term for the actions. 

It is also of interest to compare these results with the exact results for the Sauter electric field. Nikishov 0] has 
given the expected number of pairs per mode, but one must integrate over the modes to get the pair-production rate 
per area. Here we shall restrict ourselves to spinor case, as it is actually simpler. Then if one defines 

Z= 2TiqEL 2 = (61) 

the pair-production rate per area for spinors (with two spin states) can be shown after some algebra to be 



—, (qEL) 3 f cosh Zy — cosh Zx , , . , 

^ = 2 W J cosh Z- cosh Zx & ~ X )d ^ (62) 



where in this integral we have the restrictions 



— 1 < — y < x < y < 1, 

(i-z 2 )(i-y 2 )>(^) 2 = e 2 . (63) 
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This result is exact (to 1-loop in the quantum field theory) for any Z and and e < 1, pair-production being energetically 
impossible if e > 1. 

One can now show that if S <C \/l — e 2 , a good approximation to this integral expression for any e < 1 is 



For the particular subcase considered above, e -C 1 and 5 -C 1, so long as 6 <C e 4 one gets agreement with the 
instanton-method approximation given in Eq. (|58|l (after multiplying that expression by the 2 spin states of the 
spinor) . 

In particular, for S 2 <C e 4 <C 5, (where the first subdominant term in the expansion of the exponent in i|tj4[l is 
large, but not the second subdominant term), the approximation l|64(l to the exact answer (|62|l agrees with the 0-loop 
instanton-method approximation (|58|l but not with the uniform-field approximation there, which is larger by the 
factor e e ^ 4S K However, for S -C e 4 , then even the second subdominant term in the expansion of the exponent in 
(16411 is large in comparison with unity, so then the exact result is significantly different from both the uniform-field 
approximation and the 0-loop instanton-method approximation (|58|l . at least if (|64f> is a good approximation to the 
exact result to all orders in e when 5 <C \J\ — t 1 . However, we do obtain precisely l)64|l by properly calculating the 
0-loop instanton action and making a gaussian approximation for the integral over lo and . 

In the case in which e 2 « i < 1, then both approximations agree closely with the exact result. When we retain the 
assumption that £«1 but go to the opposite limit of high maximum electric field strength, so S ^> 1, then the exact 
integral expression (|62[1 reduces to the uniform-field approximation of Eq. (|60|l (after that is multiplied by the factor of 
2), rather than to the 0-loop instanton-method approximation also given there. That is, the 0-loop instanton-method 
approximation gives a result that is too large by a factor of 3ir/8 « 1.178. This is what we might expect, since the 
potential is changing rapidly so the 0-loop instanton approximation, using just the lowest-order (classical) action as 
we have done here, would not necessarily be expected to be good in this limit. And yet it is encouraging that it is 
only off by less than 20% for the pair-production rate. 

While this paper was bein g re vised, we have become aware of some related papers that have recently appeared in 
the literature [221 123. l24l 12a 12a . 1271 l28l] which use the worldline method. For a static electric field in the z-direction 
that depends only on z, the single- worldline instanton action Sq is the same as the minimum value of our 0-loop 
or classical instanton action 2S^(uj) as a function of the transverse momentum k^ and the frequency to. (The 
minimum is always at k^ = 0, and in the special cases in which scalar potential is an antisymmetric function of z, 
Aq(z) = — Aq(— z), as in, e.g., the Sauter potential, the minimum is also at u> = 0.) 

Dunne, Gies, Schubert, and Wang have now worked out the subleading prefactor contribution [28( using the worldline 
instanton method, which appears to agree with the gaussian approximation for the integrals over oj and kj_ in our 
approach using 2S^(u>). For the Sauter potential, both instanton methods give exactly the same expression as Eq. 
(I64|l which we obtained as an approximation to the exact double-integral Eq. I|t)2[) from Nikishov [lSj. 



In conclusion, using the instanton method in the Coulomb (space-dependent) gauge for constant electric and 
magnetic fields, we have obtained new forms 135|) and (|36|l for the pair-production rate in scalar QED and spinor 
QED. These formulae agree exactly with the well-known results. Finally, we suggested a generalization of these 
formulae to a constant magnetic field and an inhomogeneous electric field, and compared the results to another 
method and to the exact answer for the Sauter electric potential. 
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APPENDIX A: TUNNELING AND REFLECTION PROBABILITY 



In the context of nonrelativistic quantum theory, we derive a formula for the tunneling and reflection coefficients 
for bosons under a general potential barrier. 
Let the wave equation take the form 

y"(z) + Q(z)y(z) = 0, Q = E-V{z), (Al) 

and have three different regions: 

I = (-00, zi), Q > 0, 

II=(z l ,z 2 ), Q<0, (A2) 
JII=(z 2 , 00), Q>0. 

Here Z\ and z 2 are two turning points. We assume that an incoming wave function from —00 in region (I) is partially 
reflected back to —00 and partially tunnels through the barrier in region (II) toward 00 in region (III). Then the wave 
function in region (I) is given by 

y = y I + By* I , (A3) 

where yi has the incoming flux toward the barrier, and in region (III) 

V = Cyiu (A4) 

(assuming an outgoing flux there, which actually corresponds to an incoming group velocity for the relativistic waves 
of our paper). Here the incoming flux is normalized to unity. The flux conservation for equal momenta at z = ±00 
leads to the relation 

|B| 2 + |q 2 = l. (A5) 

Then we may find a wave function of the form [2(j 



y/ = exp 

fc=0 



>[E 5(2fc) ]> (A6) 

fc=0 

with lm(dSo/dz) > to conform with the direction of the incoming flux. The leading term in region (I) is given by 



S i0) (z) = il ^/Qizjdz, (A7) 

and in region (II) 

S {0) =i I \/~Q(z)dz (A8) 



is what may be called the 0-loop or classical action for the nonrelativistic problem. 

We may analytically continue yi from region (I) to region (III) through region (II) along a semicircle starting from 
z\ and ending at z 2 in the upper z-plane. Then z — > e~ l7r z = —z, and the leading term, for instance, takes the form 
in region (III) 



S^=-i + 1 \fQ{zj + i[ \fQ{z)+ \fW)\- (A9) 



In Eq. (|A9fl . the first term in the right hand side has the incoming flux, hence is the leading term of j/777, and the 
second term is the instanton action, a real quantity, through the barrier, and the last term gives a phase factor. We 
may thus write the analytically continued wave function as 



Vi -» e s < y if III + ay III ). (AlO) 
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So the wave function in region (I) is analytically continued in region (III) to 

VI + By* -> e s (a + B)y HI + e s (l + Ba*)y* m . (All) 
As there is no incoming wave function in region (III), we have 

e s (a + B)=C, l + a*B = 0. (A12) 
Solving Eq. I|A12() together with Eq. (|A5|) . we finally obtain the coefficients for tunneling and reflection 

\ B \ 2 = 7X^25' l c l 2 = T^s< ( A13 ) 
1 ~\~ c 1 + e Aa 

which confirms the result in Refs. @ and [IJ. 
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